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Abstract. We present the first D/H measurement in a 
damped Lya (hereafter DLA) system at z a b s = 3.025 
towards QSO 0347-3819 obtained from the UVES-VLT 
spectra. The DLA absorber has a metallicity of [Zn/H] 
w —1.25 and a relatively simple velocity structure, with 
two dominating components detected in several metal 
lines. The hydrogen Lyman series can be followed down 
to Lyl2 thanks to the high UV-Blue efficiency of UVES. 
The best fit of the Lyman series lines, and in particular of 
Lye, Ly8, LylO and Lyl2, relatively free of local contam- 
ination, is obtained when the D I absorption is included 
in the two main components. The measured deuterium 
column density yields § = (2.24 ± 0.67) x 10~ 5 close to 
other low D/H values from Lyman limit systems. The cor- 
responding values for the baryon to photon ratio and the 
baryon density derived from D/H are ij « 6 x 10~ 10 and 
ribh 2 — 0.023 respectively. 
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1. Introduction 

Deuterium (D or 2 H) is the only element entirely produced 
by nuclear reactions in the first minutes after Big Bang 
(Wagoner et al. 1967). The D yields are the most sensi- 
tive to the nuclear density at the nucleosynthesis epoch 
among the primordial light elements 3 He, 4 He and 7 Li, 
thus making the D abundance the most sensitive mea- 
surement of the baryon density in the universe (Wagoner 
1973, Schramm & Turner 1998). 

Deuterium is currently measured in the local interstel- 
lar medium (ISM), (D/H)i SM = (1.6±0.1) x 10~ 5 (Linsky 
et al. 1993), but since whenever it is cycled through stars it 
is completely burned away, extrapolation to the primordial 

* Based on public data released from UVES Commissioning 
at the VLT/Kueyen telescope, ESO, Paranal, Chile. 
Correspondence to: sdodoric@eso.org, mdessaug@eso.org, 
molaro@ts.astro.it 



D/H value requires a modeling of the Galactic chemical 
evolution. Direct D measurements of primordial material 
are thus of high interest. Adams (1976) suggested that al- 
most primordial D could be measured in low metallicity 
absorption line systems in the spectra of distant quasars 
(QSOs). This was recently achieved for a few systems, 
but with conflicting results differing by almost an order of 
magnitude. 

A few systems provide high D/H values (see e.g. Webb 
et al. (1997) who measure « 2 x 10 -4 ), while two other sys- 
tems give a low abundance at (3.39 ±0.25) x 10~ 5 (Buries 
& Tytler 1998a, 1998b). An even lower D/H estimation 
was obtained by Molaro et al. (1999), further discussed 
by Levshakov et al. (2000). Kirkman et al. (2000) mea- 
sured an upper limit of 6.76 x 10~ 5 . The handful of D 
detections found so far does not allow a firm conclusion. 
Different arguments favour a low primordial D/H ratio: 
the possible H I contamination of the D I absorption lines 
and, on the modeling side, the results by Tosi et al. (1998) 
which predict for a variety of chemical evolution scenarios 
and to be consistent with the Galactic data a maximum 
decrease of the primordial D abundance by a factor of 3. 

The paucity of suitable absorption systems for accu- 
rate D/H measurements is due to the fact that only ab- 
sorption line systems with simple velocity structures and 
with intermediate H I column densities allow the detec- 
tion of the D I lines. At too low column densities the D I 
lines are too weak for detection, whereas at high column 
densities the lines are normally washed out by the sat- 
uration of the H I line. We show here for the first time 
that in the latter case the deuterium signature can be 
successfully detected through the higher members of the 
Lyman series, when the target is a damped Lya system 
(log AT (Hi) > 20.35) at high redshift. This approach was 
first suggested by Khersonsky et al. (1995). 

2. Data 

The spectra of QSO 0347-3819 (V=17.3, z cm = 3.23) were 
obtained during the Commissioning of the Ultraviolet- 
Visual Echelle Spectrograph (UVES) on the VLT 8.2m 
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Fig. 1. Absorption metal line profiles plotted against ve- 
locity for the DLA system at z a b s = 3.025. The vertical 
scale goes from to 1 for each plotted transition. The zero 
velocity is fixed at z = 3.024856. The vertical lines mark 
the positions of 7 components. The solid thin curve rep- 
resents the best fit solution. The transitions O IA1302 and 
Si II A 1304 are from the HIRES-Keck spectra published by 
Prochaska & Wolfe (1999). 
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Table 1. Column densities 



Comp 




Ident 


logiV 
[cm" 2 ] 


b 

[km s" 1 ] 


2... 


3.024637 


Hi 


20.13 ± 0.08 


21.5 ± 0.4 






Di 


15.48 ± 0.08 


14.1 ± 0.5 


3... 


3.024856 


Hi 


20.43 ± 0.10 


23.0 ± 2.0 






Di 


15.78 ± 0.11 


16.2 ± 3.0 


7... 


3.025659 


Hi 


19.35 ± 0.05 


14.7 ± 0.4 



presented in a future paper. Here we focus on the deu- 
terium detection and the D/H ratio measurement. 

3. Analysis and results 

The absorption profiles of the DLA system at z a bs = 3.025 
are characterized by two dominating components (2 and 
3) separated by about 20 km s _1 with the red one slightly 
stronger than the blue as it can be seen from the non- 
saturated metal lines in Fig. [l| The strong and satu- 
rated metal lines reveal that additional material in smaller 
amount is present redwards the two main features (com- 
ponents 4, 5, 6 and 7). The region on the blue side of the 
main absorption components is sharp and relatively free 
from material with only one weak component (component 
1) at about —30 km s . In total seven components are 
needed to fit the metal lines absorption profiles, with two 
main components containing about 80% of the total col- 
umn density per transition. 

The DLA system at z a b s = 3.025 is a very good can- 
didate for the deuterium analysis since it shows a rela- 
tively simple velocity structure dominated by two strong 
components, a low metallicity of [Zn/H]Q « —1.25, indi- 
cating that the measured D/H will be representative of 

1 [X/H] = log(X/H) -log(X/H) 



Kueyen telescope at Paranal, in 1999. The instrument is 
described in D'Odorico et al. (2000). Two exposures of 
4500 s each, covering the spectral range from 3650 to 4900 
A and from 6700 to 10000 A were obtained with a resolu- 
tion of « 6.9 km s _1 and ps 5.7 km s _1 respectively. The 
individual spectra were reduced using the UVES data re- 
duction pipeline implemented in the ESO MIDAS package. 
The final spectrum reaches a S/N varying from 20 to 35. 

QSO 0347-3819 shows a damped Lya system at z a b s = 
3.025 which has been studied in detail by Centurion et al. 
(1998), Ledoux et al. (1998) and from HIRES-Keck obser- 
vations by Prochaska & Wolfe (1999). The UVES obser- 
vations however provide the first high-quality data in the 
UV (A < 4900 A) and in the near-IR (A > 6700 A). They 
allow the abundance measurements of new features such 
as O, P, Ar and Zn, in addition to the N, S, Si and Fe 
abundances measured in the previous studies. Full abun- 
dance analysis of the DLA system at z a h s = 3.025 will be 
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Fig. 2. Damped Ly/3 absorption profile. The vertical lines 
mark the positions of the three major hydrogen compo- 
nents (see text for details) . The solid thin curve represents 
the best fit solution and the dotted thin curves delimit the 
range of uncertainty of the hydrogen column density. 
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Fig. 3. Lye, Ly8, LylO and Lyl2 profiles plotted against 
velocity. The vertical scale goes from to 1 for each plotted 
transition. The zero velocity is fixed at z = 3.024856. The 
vertical lines mark the positions of the three hydrogen 
components. The solid thin curve represents the best fit 
solution. 

primordial D/H and the hydrogen Lyman series lines can 
be followed down to Lyl2. 

We used a x 2 minimization routine f itlyman (Fontana 
& Ballester 1995) in MIDAS to fit Voigt profiles to the ob- 
served absorption profiles, and obtain for each fitted ab- 
sorption component the wavelength, the column density 
N, the Doppler parameter b and the corresponding errors. 

In the case of high hydrogen column densities like in 
DLA systems, we expect the neutral and low ionization 
metal lines to trace the H, therefore we model the Lyman 
series (Ly/3 to Lyl2) absorption profiles with three fea- 
tures corresponding to the two main components (2 and 
3) and the reddest component (7) as determined from the 
metal lines (Fig. Q). The contribution of the weaker com- 
ponent 7 is required only to better constrain the fit on the 
red edge of the Lyman lines. The relative intensities of 
these three major hydrogen components were scaled with 
the metal lines assuming they have approxirnatively the 
same abundance ratios from component to component. 

Starting from this basis we obtained the final H I col- 
umn densities and 6-values (free parameters) by fitting 
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Fig. 4. Same as Fig. [| The fit in this case is performed 
by considering the two main components for D and the 
three components for H (see text for details). 

simultaneously the lines of the Lyman series and by as- 
suming the same redshift for the three H components as 
for the metal lines. The Ly/3 absorption profile provides a 
very good constraint to the total H I column density and 
the Ly8, Ly9, LylO and Lyl2 profiles (which are the ones 
free from strong contamination) to the ^-values. From the 
best fit (Fig. §, Fig. |, Table |) we derived a total Hi 
column density of log iV(Hi) = 20.63 ±0.09 in close agree- 
ment with the Pettini et al. (1994) value of 20.7 ± 0.1. 

Adding other components to the hydrogen fit with 
low column densities as the ones observed in the stronger 
metal lines, component 1 at —30 km s _1 and components 
4, 5 and 6 at slightly higher redshifts than the two main 
components, does not change significantly neither the to- 
tal H I column density of the two main components nor the 
fit on the blue and red wings of the Lyman lines. Fig. |^ 
however clearly shows that the fit with three hydrogen 
components systematically fails to reproduce the edge of 
the blue wing in the higher members of the hydrogen Ly- 
man series, Lye, Ly8, LylO and Lyl2, at about —82 km 
s _1 from the two H main components (2 and 3) which is 
the expected displacement of the corresponding D lines. 

We then add to the model the deuterium by assuming 
its contribution only to the two main components (2 and 



4 



S. D'Odorico, M. Dessauges-Zavadsky and P. Molaro: Deuterium in a damped Lya system 



3). In the fitting procedure we assumed the same redshift 
for the two fitted components as for H and we left the 
column densities and the broadening parameters b as free 
parameters. 

The best fit performed over D I Ly8 and LylO improves 
the x 2 °f th e hydrogen fit by a factor of 3. It reproduces 
nicely the blue wing of the Ly8, LylO and Lyl2 absorption 
profiles and is consistent with Lye (see Fig. [|). It gives a 
deuterium abundance of 



D 
H 



= (2.24 ±0.67) x 10" 



The error follows from the H and D fits and it includes 
the uncertainties in placing the continuum level and the 
errors on the b- values. The D I column densities and the 
b- values of the two fitted components are given in Table § 
To check the stability of our result when the assump- 
tion that H and D mimic the metal structure is relaxed, we 
made a fit of the Ly/3, Lye, Ly8, LylO and Lyl2 lines with 
three H and D components by keeping both their redshifts 
and column densities as free parameters. We obtained a 
satisfactory fit with shifts of the components with respect 
to the values given in Table [l] of less than 4 km s _1 . The 
derived D/H ratio was comprised within our D/H error 
range. A fit with only one strong H component and the 
corresponding D component gave also a D /H ratio within 
our D/H error range. 

4. Discussion 

Three arguments support the interpretation of the ab- 
sorption feature detected on the blue wing of some of the 
higher members of the Hi Lyman series, Lye, Ly8, LylO 
and Lyl2, as the Di absorption with the same two main 
components as detected in the metal lines of the DLA 
system. 

First, the absorption in the blue wing of the Lyman se- 
ries lines could also be explained by a hydrogen interloper 
associated with the damped Lya system with a column 
density between 10 15 7 and 10 16 2 , a 6-value between 15 
and 25 km s _1 and placed between —80 and —100 km 
s _1 . However, from the density distribution of the Lya 
clouds in the forest at z ~ 3 (see e.g. Kim et al. 2000), 
the probability to have such a cloud at that position is 
smaller than 1/1000. The lack of any metal component 
at this velocity in the strong and saturated lines of the 
DLA system (see Fig. [l]) provides an additional evidence 
for discarding this possibility. 

Secondly, the contamination on the blue wing of the 
Lyman series lines by different H I interlopers at different 
redshifts which would mimic the same deuterium abun- 
dance as derived from the Di Lye, Ly8, LylO and Lyl2 
lines is even more unlikely. 

Finally, the D and H column densities and broadening 
parameters b resulting from the fits are consistent one rel- 
atively to other in the two fitted components (2 and 3): 



the derived D/H ratios are the same and the b(D i)/6(H i) 
ratios are close to what is expected in the thermally dom- 
inated case. 

This measurement of D is the first made in a DLA 
system. It shows that the DLA systems with their low 
metallicity ISM are a very promising class of absorbers for 
measurements of the D/H ratios at high redshifts, when 
it is possible to measure the higher members of the Ly- 
man series. A systematic program of measurements using 
UVES data is under way. 

The derived D/H ratio of (2.24 ± 0.67) x 10~ 5 is close 
to the low values obtained by Buries & Tytler in Lyman 
limit systems (1998a, 1998b) and it makes the claim of 
the primordial low D/H ratio more robust. Taken at face 
this ratio gives a baryon to photon ratio, rj = nb/rij, of 
w 6.3 x 10~ 10 (Buries et al. 2000). This -q implies a he- 
lium abundance (in mass) of Y p = 0.2480 and a lithium 
abundance of ~ 5 x 10~ 10 which are both larger than 
presently allowed by observations of He in extragalactic 
Hn regions and of Li in halo stars (Izotov & Thuan 1998, 
Bonifacio & Molaro 1997). On the other hand an rj of 
w 6.3 x 10~ 10 corresponds to a present-day baryon density 
of n b h 2 = 0.023FI which remains significantly lower than 

1 0.009 
0.008 

(Jaffc et al. 2000). 



the fl b h 2 = 0.032t°;^ derived from CMB anisotropy 
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